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Abstract 

First, we construct the Taub-NUT/bolt solutions of (2fc-|-2)-diniensinaI Einstein-Maxwell gravity, 

^S| ! when all the factor spaces of 2A;-dimensional base space B have positive curvature. These solutions 

> 

depend on two extra parameters, other than the mass and the NUT charge. These are electric 



o 
o 



charge q and electric potential at infinity V . We investigate the existence of Taub-NUT solutions 

and find that in addition to the two conditions of uncharged NUT solutions, there exist two extra 

conditions. These two extra conditions come from the regularity of vector potential at r = A^ and 

the fact that the horizon at r = A^ should be the outer horizon of the NUT charged black hole. We 

^ . find that the NUT solutions in 2A;-|-2 dimensions have no curvature singularity at r = A^, when the 

^_ ■ 

2A;-dimensional base space is chosen to be CP . For bolt solutions, there exists an upper limit for 

K> ' the NUT parameter which decreases as the potential parameter increases. Second, we study the 

thermodynamics of these spacetimes. We compute temperature, entropy, charge, electric potential, 

action and mass of the black hole solutions, and find that these quantities satisfy the first law of 

thermodynamics. We perform a stability analysis by computing the heat capacity, and show that 

the NUT solutions are not thermally stable for even fc's, while there exists a stable phase for odd 

fe's, which becomes increasingly narrow with increasing dimensionality and wide with increasing V . 

We also study the phase behavior of the 4 and 6 dimensional bolt solutions in canonical ensemble 

and find that these solutions have a stable phase, which becomes smaller as V increases. 
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I. INTRODUCTION 

It has been known for quite long time that the area law of entropy holds for black 
holes or black branes in any dimension 2k + 2, when the U{1) isometry group associated 
with the (Euclidean) timelike Killing vector dr has a fixed point set on the surfaces of 2k 
codimension [l|. However, in recent years, it has been shown that entropy can be associated 

nn 

with a more general class of spacetimes [2;|3|]. In these spacetimes, the f/(l) isometry group 
can have fixed points on surfaces of any even codimension, and the spacetime need not be 
asymptotically fiat or asymptotically anti-de Sitter (AdS). In this more general class, the 
entropy is not just a quarter the area of the 2/i;-dimensional fixed point set. Such situations 
occur in spacetimes containing NUT-charges. In 2A; + 2 dimensions they not only can have 
2fc-dimensional fixed point sets (called bolts), they also have fixed point sets with dimension 
less than 2k (called NUT). Here the orbits of the U{1) isometry group develop singularities, 
which are the gravitational analogues of Dirac string singularities and are referred to as 
Misner strings. When the NUT charge is nonzero, the entropy of a given spacetime includes 
not only the entropies of the 2/i;-dimensional bolts, but also those of the Misner strings. In 
a very basic sense, gravitational entropy can be regarded as arising from the Gibbs-Duhem 
relation applied to the path- integral formulation of quantum gravity '4']. If one calculates 
the finite total action / evaluated on the classical solution, then the entropy may be written 
as 

S = f3{M - r,Q - I (1) 

where Cj and Fj are the conserved charges and their associate chemical potentials respectively. 
The original asymptotically locally flat NUT/bolt solutions in four dimensions have been 
constructed in Ref. j5j. There are known extensions of the Taub-NUT/bolt solutions to 
the case when a cosmo logical constant is present. In this case the asymptotic structure is 
only locally de Sitter (for positive cosmological constant) or AdS (for negative cosmological 
constant) and the solutions are referred to as Taub-NUT-(A)dS metrics. Generalizations to 



higher dimensions follow closely the four-dimensional case |6l.l7l.l8Ll9LllO. 
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. Taub-NUT 



solution of the Einstein equations with multiple NUT parameters has been investigated in 



[13[. Also, charged Taub-NUT solution of the Einstein-Maxwell equations in four dimensions 
is known [1^ , and its generalization to six dimensions has been done in Ref. |l5|, |l6| . The 
existence of NUT charged solutions of Einstein- Yang-Mills and Einstein- Yang-Mills-Higgs 



theory and their thermodynamics have also been considered [17|. Recently, the existence 
of Taub-NUT/bolt solution in Gauss-Bonnet and Gauss-Bonnet-Maxwell gravity have been 
studied by one of us jl8|, |l9[. In this paper we, first, construct the {2k + 2) -dimensional 
NUT/bolt solutions of Einstein gravity in the presence of electromagnetic field and, second, 
we calculate their conserved charges through the use of the counterterm method and analyze 
the thermodynamic behavior of these electrically charged NUT/bolt solutions. 

The outline of our paper is as follows. We give a brief review of the counterterm method 
in Sec. |n] In. Sec. IIIII we construct the {n + l)-dimensional Taub-NUT/bolt solutions of 
the Einstein-Maxwell gravity. The thermodynamics of electrically charged NUT solutions 
in 4-dimensions is investigated in Sec. II VI and that of 6-dimensional solutions is considered 
in Sec. [V] Generalization of these subjects to the case of {2k + 2)-dimensional solutions is 
done in Sec. IVII We finish our paper with some concluding remarks. 



II. GENERAL FORMALISM 

The gravitational action for Einstein gravity in {n + 1) dimensions in the presence of 
cosmological constant and electromagnetic field is 

^o = -T^ [ d^^'x,/^{n - 2A - F'^^F^,) -^ I rf"xv/l7^(7), (2) 

where A = —n{n — l)/(2/^) is the cosmological constant, TZ is the Ricci scalar, F^j, = 
dfj^Aij — d^A^ is electromagnetic tensor field and A^ is the vector potential. The first term 
is the Einstein-Hilbert action and the second term is the Gibbons Hawking boundary term 
which is chosen such that the variational principle is well-defined. The manifold Ai has 
metric g^i, and covariant derivative V^. K is the trace of the extrinsic curvature K'^'^ of any 
boundary(ies) dAi of the manifold Ai, with induced metric(s) 7jj. 

In order to obtain the Einstein-Maxwell equations by the variation of the volume integral 
with respect to the fields, one should impose the boundary condition 6A^ = on dAi. 
Thus the action Q is appropriate to study the grand-canonical ensemble with fixed electric 
potential J20(]. To study the canonical ensemble with fixed electric charge one should impose 



the boundary condition 6{n°'Fab) = 0, and therefore the gravitational action is 21 1 
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^g = Ig-^ I (Tx^-^UaF'^'A,, (3) 



where n^ is the normal to the boundary dA4. Varying the action Q or Q with respect 
to the metric tensor g^^iy and electromagnetic tensor field F^^,, with appropriate boundary 
condition, the equations of gravitation and electromagnetic fields are obtained as 

G'..-^^^^^.. = 87rT(r), (4) 

V^F'^^ = 0, (5) 

where G^i, is the Einstein tensor and Tf^u = 2F'' Fpy — ^Fp„FP"g^y is the energy-momentum 
tensor of electromagnetic field. 

In general the action Iq or Iq are diverged when evaluated on solutions, as is the Hamil- 
tonian and other associated conserved charges. One way of eliminating these divergences 



is through the use of background subtraction [22, 12^, in which the boundary surface is 
embedded in another (background) spacetime, and all quasilocal quantities are computed 
with respect to this background, incorporated into the theory by adding to the action the 
extrinsic curvature of the embedded surface. Such a procedure causes the resulting physical 
quantities to depend on the choice of reference background; furthermore, it is not possible 
in general to embed the boundary surface into a background spacetime. For asymptotically 
AdS solutions, one can instead deal with these divergences via the counterterm method 
inspired by AdS/CFT correspondence J2J|. This conjecture, which relates the low energy 
limit of string theory in asymptotically AdS spacetime and the quantum field theory on its 
boundary, has attracted a great deal of attention in recent years. The equivalence between 
the two formulations means that, at least in principle, one can obtain complete information 
on one side of the duality by performing computation on the other side. A dictionary trans- 
lating between different quantities in the bulk gravity theory and their counterparts on the 
boundary has emerged, including the partition functions of both theories. In the present 
context this conjecture furnishes a means for calculating the action and conserved quantities 
intrinsically without reliance on any reference spacetime |25l l26l |27 | by adding additional 
terms on the boundary that are curvature invariants of the induced metric. Although there 
may exist a very large number of possible invariants one could add in a given dimension, 
only a finite number of them are non vanishing as the boundary is taken to infinity. Its many 
applications include computations of conserved quantities for black holes with rotation, vari- 
ous topologies, rotating black strings with zero curvature horizons and rotating higher genus 
black branes J28| . Although the asymptotic structure of the NUT charged solutions of Ein- 



stein gravity in the presence of a cosmological constant is only locally (A)dS, the counterterm 
method has been employed to calculate the action and the conserved charges of them. The 
conserved quantities of asymptotically locally AdS solutions in four dimensions have been 
calculated in [29|, those of higher- dimensional solutions in J30|, and those of asymptotically 
locallv dS solufous in Q. Also theconserved quantises of the asymptotica.lv .ooally AdS 
solutions have been calculated in |32]. Usually, the counterterm method applies for the 
case of a specially infinite boundary, but it was also employed for the computation of the 
conserved and thermodynamic quantities in the case of a finite boundary 33|. The general 
correspondence formula is [2J| 

f I?^e-^-4'*^-[*l = (exp I (rxO{x)^o{x)), (6) 

where the functional integral on the left hand side is over all the fields \l/ whose asymptotic 
boundary values are \l/o, and O denotes the conformal operators of the boundary conformal 
field theory. In the classical limit, the correspondence formula can be written as [341 ? ] 

lAds[^o] = WcFxl^oi (7) 

where lAds is the classical on-shell action of an AdS field theory, expressed in terms of 
the field boundary values \l/o, and Wcft is the CFT effective action. However, one should 
expect I AdS to be divergent as it stands, because of the divergence of the AdS metric on 
the AdS horizon. Thus, in order to extract the physically relevant information, the on-shell 
action has to be renormalized by adding counterterms, which cancel the infinities of Iq in 
the absence of matter. This counterterms up to nine dimensions is J2^ 

1 /■ .n.. ^,n-l iein-3) 



-rtafc-K — TT TT-n^ + 7 7:TT7 ITT ;T 77 TTJ^J^abJ^ 



2{n-A){n-2y\ 4(n - 1) J (n - 2)3(n - 4)(n - 6) '4(n - 1) 

rifr, J- 9] n — 1 1 

^ ^ ^ R^-2R^'R,mR'"' + -, -R''VaV,R-R^'UR,, + - -RUR]} {^) 



16(n-l)2 ""'•""" ' 2(n-l) " " "" ' 2(n - 1 

where R, Rated, and Rat are the Ricci scalar, Riemann and Ricci tensors of the boundary 
metric jab, and Q{x) is the step function which is equal to one for x > and zero otherwise. 
In the presence of matter, one may encounters with some divergencies. For n > 4 the 
electromagnetic field will cause a power law divergence in the action which can be removed 



by a counterterm of the form 
I^ = ±fr.^ |e(n-5)(n-8) Pe(n-7) (5n - 11) 



SttJ ^ ' [ 32 (n-4) 4 '48(n - l)2(n - 2)(n - 6)' 

I (7n-66) pfec I (^-8) /^ pa6N2 

+ 12(n-6)(n-2) " "'^ +12(n-4)2^^"^ ^ 
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-F'^^(VfeV'=F,,-VaV^Fefe)]| (9) 



A7{n-A){n-Qy 

Thus, the total finite action can be written as a hnear combination of the gravity term (j21) 
and the counterterms (jH)) and Q- Having the total finite action, one can use the Brown and 
York definition of energy-momentum tensor J22] to construct a divergence free stress-energy 
tensor. This tensor is 

1 n — 1 1 1 

Svr / n — 2 2 

/39(n-5) 1 ab.pcdp ^^^^r2x ^ Dpgfe 

+ 7 TV? o\2 ~o^ '--^ ^<"^'~T( TT-^ '' ~ To n\^^ 

[n — 4:)[n — ly 1 4(r2 — 1) \ln — 1) 

n — 1 1 r^P™ 

+2R,dR''''"' - -, --V"V'i? + V'i?"' - 7T7"'V'i?] + ...} + -^ (10) 

2(n - 1) 2[n - 1) d^iab 

The explicit form of the stress-energy tensor due to the electromagnetic field counterterm 
will be given in Sec. |3 To compute the conserved charges of the spacetime, we choose a 
spacelike surface S in dM. with metric (Tij, and write the boundary metric in ADM form: 

'yatdx^dx" = -Af^dt^ + Oi^ {d^' + V'dt) [d^^ + VUt) , (11) 

where the coordinates </?* are the angular variables parameterizing the hypersurface of con- 
stant r around the origin, and M and V* are the lapse and shift functions respectively. The 
conserved charges associated to a Killing vector E,'^ is 

Q(0 = j d^-'xy/^u'^Tate, (12) 

where a is the determinant of the metric aij and m" is the normal to the quasilocal bound- 
ary hypersurface S. For boundaries with timelike Killing vector {C, = dt) one obtains the 
conserved mass of the system enclosed by the boundary S. In the context of AdS/CFT 
correspondence, the limit in which the boundary E becomes infinite (Sqo) is taken, and the 
counterterm prescription ensures that the action and conserved charges are finite. No em- 
bedding of the surface E in to a reference of spacetime is required and the quantities which 
are computed are intrinsic to the spacetimes. 



III. {2k + 2)-DIMENSIONAL TAUB-NUT/BOLT SOLUTIONS IN EINSTEIN- 
MAXWELL GRAVITY 

The Euclidean section of the {2k + 2) -dimensional charged Taub-NUT/bolt spacetime 
can be written as 

ds^ = F{r) {(It + NAf + F-\r)dr^ + (r^ - N^)dEB, (13) 

where r is the coordinate of the fibers S^ and A is the Kahler form of the base space B, 
N is the NUT charge and F{r) is a function of r. The metric dEjs is a 2/c-dimensional base 
space Einstein-Kahler manifold B. 

Here, we consider only the cases where all the factor spaces of B have positive curvature. 
Thus, the base space B may be the product of CP spaces for all values of k. For com- 
pleteness, we give the 1-forms and the metrics of these factor spaces. The 1-forms and the 
metrics of CP is 

Ak = 2{k + 1) sin" a(#fc + ^^fc-i)' (14) 



dT^i = 2{k + l)Uek + ^in'ikCos'Udi^k + ^Au-iY + —sin'ikdT.i_A (15) 

where Ak-i is the Kahler potential of CP ~^. In Eqs. p4|) and p5|) ^k and "ipk are the extra 
coordinates corresponding to CP with respect to CP ^^. The metric CP is normalized 
such that, Ricci tensor is equal to the metric, R^^, = g^^^. The 1-form and the metric of CP^ 

{S^) is 

Ai = 4sin2^i#i (16) 

ciSi^ = 4(rfei^ + sin2^icos^6^^i^) (17) 

and those of CP''' can be constructed through the use of Eqs. (IT^ and (IT^ . 

A = h{r){dT + NA), (18) 

where h{r) is a function of r. The electromagnetic field equation (jSJ for the metric ()1H|1 with 
vector potential (fTHj) is 

(r^ - N^fh"{r) + 2kr{r^ - N^)h'{r) - AkN^h{r) = (19) 
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where prime denotes a derivative with respect to r. The solution of Eq. ()19|) may be 
expressed in terms of hypergeometric function 2-^1 ([a, b], [c], z) in a compact form as 

Kr) = pr^ + ^(2^ - 1)^" 2^i([-^' -^]' t^]' ^) (20) 

where ^ and q are two arbitrary constants which correspond to electric potential at infinity 
and charge respectively. To find the function F{r), one may use any components of Eq. (^. 
The simplest equation is the tt component of these equations which can be written as 



(21) 



with the solution 



+k(2k + l)(s^ - NY + kl^^' - N^ ds - j^TZJpf. (22) 

where h{r) is given in Eq. ((201) • One may note that the above solution reduces to the 
(2/i; + 2)-dimensional solution given in [30| in the absence of electromagnetic field (h{r) = 0). 



A. NUT Solutions 

The solutions of Eq. ^U^ describe NUT solutions, if 

(I) F(r+ = A^) = 0. 

(II) /3 = An/F'{r+) = 47r(A; + 1)N 

(III) h{r+ = N) = 0. 

(IV) F{r) should have no positive roots at r > A^. 

The first condition comes from the fact that all the extra dimensions should collapse to zero 
at the fixed point set of d/dr, the second one ensures that there is no conical singul arity 
with a smoothly closed fiber at r = A^. Of course, if one use the Misner's argument J37|, 
the period of time coordinate is found to be different from Anik + 1)A^ when the space is 

pn 

singular [7|, |15|| . Here we consider the period of Euclidean time coordinate by the elimination 
of conic singularities. The third condition comes from the regularity of vector potential at 
r = N and the fourth one comes up since r = N should be the outer horizon. Condition 
(III) gives a relation between q and V as 



Using the first two conditions with Eq. ()23p . one finds that Einstein- Maxwell gravity in even 
dimensions admits NUT solutions with any base space when the mass parameter is fixed to 
be 

m„ = N^''-'{2B{k)[P - 2{k + 1)N^] + 4:D{k)V^} (24) 

where Bk and D{k) are 

r(| - A;)r(A; + 1) ^ , 

Finally, we should apply the fourth condition. As we will see in the next section, the fourth 
condition is not necessary in 4 dimensions. However, this condition restrict the value of Vn 
to be less than a critical value Vcrit in dimensions higher than four. To find Vcrit we proceeds 
as follows. We define the function gnntir) as the numerator of -Fnut(^) = F{q = g„,m = 
m„, r)/{r — N) which is positive at r = A^, and solve the system of two equations 

9Ur) = ^27) 

for the unknown V and r. The V obtained by this method is the critical value V^rit- We 
consider the solutions of system of two equations P7j) in the following sections for various 
dimensions. 

Now, we calculate the electric charge of the solutions. To determine the electric field we 
should consider the projections of the electromagnetic field tensor on special hypersurfaces. 
The normal to such hypersurfaces is 

u' = ^, u^ = 0,u^ = ^, (28) 

where A/" and V* are the lapse and shift functions respectively. The electric field is E'^ = 
g^'^FpyU'^ , and the electric charge can be found by calculating the fiux of the electric field at 
infinity, yielding 

Qk = {2k-l){ATxf-\. (29) 

The electric potential $, measured at infinity with respect to the horizon, is defined by J3^ 

$ = ^mX^ Ir^oo - ^mX1,=,^ , (30) 

where x = f^r is the null generator of the horizon. We find 

$„ = (-1)'+V„. (31) 



B. Bolt Solutions 

The conditions for having a regular bolt solution are (I) F{r = Vb) = 0, (II) F'{rb) = 
[{k + 1)A^]~-'^ and (III) h{ri,) = with r^, > A^. Condition (II) follows from the fact that we 
want to avoid a conical singularity at the bolt, together with the fact that the period of r 
will still be 47r(/c + 1)A^. The first and third conditions gives 

'^'-■^^[^) (2^-2.-1) + P h\ ' ) {2k-2^ + l) 

+ / SL ILI [(^2 _ jv2)2/,'2(^) _ AN^kh^ir)]dr. (32) 

J kr^ 

where 14, is the solution of h{rh) = 0. The second condition together with Eq. (jH^ gives the 

following equation for r?, 

kfrl - k{k + l)N[{2k + l){rl - N^) + f]rl - {k + l){2k - iyfV,^N{rl - N^) = 0. (33) 

Equation (J33|l for fixed values of / and VJ, is a quartic equation in r?, and has two real solutions 
greater than A^ provided N < A^max and one real solution greater than r^ if A^ = A^max- 
Numerical calculation shows that A^max decreases as V increases in various dimensions. We 
will give the equation which is satisfied by A^max in four and six dimensions in Sees. IIVI and 
E respectively. The electric charge is the same as the NUT solutions, Eq. ()29p. while the 
electric potential may be obtained by use of Eq. (jHUj) as 

$. = tl)V^ ^ (34) 

{2k-l)N'>' ,F,{[-l-k],[l],^) 

IV. THERMODYNAMICS OF 4-DIMENSIONAL SOLUTIONS 

In four dimensions, the functions h{r) and F[r) given in Eqs. (j2m) and (j22j) are 

"(O = ^r^ + f S- (35) 

nr) = /2(,. _ ;v2) (36) 

A. NUT Solutions 

Since m„ and g„ of Eqs. (j^ and (PH|) become 

mn = 2Nl~\P-AN^), (37) 
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qn = -2iVK, (38) 

the function F{r) for NUT solution can be written as 

r'^{r + N)+ P{r - N) - N'^{5r - 3N) 
"^"■^ ~ (r + N)P • 

It is notable that Fn{r) is independent of q and V, and therefore there is no restriction on 
Vn from the fourth condition. The electric charge and potential are 

Q = q, (39) 

^n = -— (40) 

Using Eqs. and (jHI), the Euclidean actions in the grand-canonical and canonical 
ensembles can be calculated as 

J„ = ^(z2_2Ar2 + 2/V2), (41) 

4 = ^^(/' - 2N^ - 2/V2) (42) 

Also the total mass may be calculated as 

Mr. = ^, (43) 



and as in the case of uncharged solutions j30|] the total angular momentum is zero. The 
Gibbs and Helmholz free energies are 

GiT ^) = L = m^' + ^r-'T^-^ (44) 

J _ 1 1 - 32/^7r^r^ + 1024g^/^7r^T^ 

where T is the Hawking temperature which is inverse of /3„ = SttA^. The entropy can be 
obtained as 

The entropy can also be obtained through the use of Gibbs-Duhem relation (Q), where Ci = Q 
and Fj = $. As one can see from Eq. PUJ) . the entropy is positive for 1) t^ > GN"^ and also 
2) /^ < GN"^ provided \Vn\ > {2iy^y/12N'^ — 2P. It is worth to mention that the second case 

11 



occurs only in the case of 4-diniensional NUT solution in the presence of electromagnetic 
field. 

The Smarr-type formula for mass versus extensive quantities S and Q may be written as 

P 
where Z is 



M„ = ''"^'r"'^ (47) 



2 = ^(1 ± VP + 12Q2 - eS/vr) (48) 

One may then regard the parameters S and Q as a set of extensive parameters for the 
mass J^{S,Q) and define the intensive parameters conjugate to S and Q, which are the 
temperature and the electric potential respectively. One obtains 

'dMn\ 1 



fdM\ q 

for + and — sign in Eq. ()48|) provided 12N'^ > P and 12N'^ < P respectively. Equations 
PUj) show that the the quantities T„ and $„ coincide with those which was calculated in 
Sec. flllljl . Thus the thermodynamic quantities calculated in Sec. pilll for Nut solution, 
satisfy the first law of thermodynamics dM = TdS + (^dQ. 

The stability of a thermodynamic system with respect to the small variations of the 
thermodynamic coordinates, is usually performed by analyzing the behavior of the entropy 
near equilibrium. The local stability in any ensemble requires that S{Y) be a concave 
function of its extensive variables or that its Legendre transformation is a convex function 
of the intensive variables. The stability can also be studied by the behavior of the energy 
M(X) which should be a convex function of its extensive variable. In our case the mass 
M (5*, Q) is a function of entropy and charge. The number of thermodynamic variables 
depends on the ensemble that is used. In the canonical ensemble, the charge is a fixed 
parameter, and therefore the positivity of the heat capacity Cq = T{dS/dT)Q is sufficient 
to ensure local stability. The heat capacity Cq at constant charge is 

Cr. = ^{12N'-f) (50) 

which is positive provided P < 12N'^. Thus, the NUT solution is stable in two cases. The 
first case is when 6N^ < P < 12N'^ for which both the entropy and the heat capacity are 
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positive. The second case is when P < GN"^ provided |Ki| > {2l)^^y/12N'^ — 2P. This second 
case is a stable phase for charged NUT solutions , which does not occur in the absence of 
electromagnetic field. That is, by turning on the electromagnetic field, the stable phase of 
the four-dimensional NUT solution gets a new zone. 

B. Bolt Solutions 

In 4 dimensions, the bolt mass, charge and electric potential may be written from Eqs. 
dSl, dSSl), and dSl as 

_ ^4+(/2_eJV2)^2 + ^2(^2_3^2) y2(^2_^2)2 

.. - -^-^^, (52) 

where Vf, is the solution of the following equation 

6Nrt - t^rl + 2N{f - 3A^^ + fV^Yl - 2fV'^N^ = 0. (54) 

For fixed values of / and V, Eq. (J5l|l has two real solutions greater than A^ provided 
A^ < A^max and one real solution rj, > N for A^ = A^max where A^max is the smallest root of 
the following equation 

62208Ari° - 82944(1 - V^)fN^ + {AU72V^ - 276481^^ + 41904)/^Ar'' 

+ (9216V"^ + 9216V"^ - AAQAV^ - 9648)/^ A^^ + (768\/*^ + 3072^*^ 

+30241/^ + 16321^2 + 912)/^ A^2 - (16^^ + 48F^ + 211^^ + 16)/^° = (55) 

Thus, we have bolt solution for A^ < A^max- For example, for / = 1 and V^ = 1, the maximum 
value of A^max is 0.1033. It is a matter of numerical calculation to show that A^max decreases 
as V increases. The finite actions in grand-canonical and canonical ensembles are 



r 



— TT- 



t'rl + Ar2(3Ar2 - /2)]^2 _ (^4 ^ 4JY2^2 _ ]V4)/V2 



(3r2 - 3A^2 + /2)^2 ^ (^2 _ ]V2)/2V2 



(56) 



\rl - l\l + N\ZN'' - /2)]r2 + (3r,^ + N^)l''V^ 

'^ {^rl - 3Ar2 + /2)^2 ^ (^2 _ ]V2)/2y2 ' ( ^) 



13 



200- 


; 












150- 


\\ 














100- 


\ \ 














50- 


\ 


\ 












0" 




""^ 


^~~~^ 


-^ 


— 


J 




0.04 


0.06 


0.08 


0.1 


0.12 


0.14 





FIG. 1: Entropy (bold line) and heat capacity (solid line) of 4-diniensional bolt solution for / = 1 
and y = in the allowed range of N. 

while the entropy and specific heat are calculated as 



S, 



TT 



(3r2 - 3Ar2 + i2y2 ^ ^^2 _ jv2)Z2y2 "^b 
(3r2 - 3Ar2 + i2y2 + (^2 _ Ar2)/2^2' 



+7r 



(58) 
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Cb = {[4:N{rl + N^)fV^ + ?,QNrt - hl^r'l + ANrl{2l^ - ?,N^) - /^iVVj/V,^}"' x 
2Nt^[- {ASN^rll^ + QNrll"^ + QNY + ASN^ll"^ - 108 A^V^/^) V^ - {264N^rtf 
+36Nrl°P + l2Nrll^ + 216A^^/V^ - AhN^rll^ + mN^t^rl - 76N^4l^ - br^ 
+5N^rll^ - IGN'^rll^ - dGONH'^rf + 45N^rll^ + 112N\^^l^)V^ + 264N\IP 
+ l5rlH^ + 360A^V^ - 213iVV^/^ - 36Nrl^l^ - 396ArV^ - QNrfl'^ + 72N'^f 
-12A^5/V^ - lON^r^l^ + brfl^ + IN^rll"^ - Q3NHV^, + 1368N^rl^ - 5ANrl 
-126N\^ - IGN^r^J^ + UN'^rll^ - 27 NH^rl]. 

In Figs. [T[ 121 and 01 the entropy and specific heat are plotted for / = 1, V = 0, \^ = 1 
and l^ = 3 up to A^max for outer horizons of the bolt solutions respectively. As one can see 
from these figures, the bolt solution is stable in the range < A^ < A^st < A^max? while A^st 
is less than A^max for ^ 7^ 0, and equal to A^max for ^ = 0. That is, the uncharged solution 
is stable in the whole allowed range < A^ < A^max, while the stable phase becomes smaller 
as V increases. 
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FIG. 2: Entropy (bold line) and heat capacity (solid line) of 4-diniensional bolt solution for / = 1 
and y = 1 in the allowed range of N . 




FIG. 3: Entropy (bold line) and heat capacity (solid line) of 4-diniensional bolt solution for / = 1 
and y = 3 in the allowed range of N . 

V. THERMODYNAMICS OF 6-DIMENSIONAL SOLUTIONS 



In six dimensions the base space B can be the 4-diniensionaI space CP^ or the product 
of two 2-dimensional spaces S"^ x S"^ . From Eqs. (|2m) and (|22p the functions /i(r) and Fir) 
may be written as 

Kr) = T-T^^-TT^. - ^ ... .r... -■ (59) 



(^2 _ p^2y 



(^2 _ p^2y 



Fir) 



3^6 ^ (;2 _ i5]v2)r4 + 3N\15N^ - 2Py + 3N\5N^ - P 

/2(r2 - my 

3r2-Ar2 ^6 ^ ^5Jy2^4 _ g^4/^2 _ ^2N 

a — AN ^ -V — 

2{r^-NY^ [r^-my 



mr 



(^2 _ JY2-)2 

IQN'^r^ 

-qV. (60) 



(r2-iV2)4 



15 



A. NUT Solutions 

The mass and charge parameters of the NUT solutions may be obtained from Eqs. 
and (pH|) as 



m„ = -::-^(r^-6N' + 9PV^), (61) 

g„ = 8iW„. (62) 



Substituting the above parameters in Eq. (jHn|l one obtains 

Fn{r) = 3(l.Vi^!^2 [-12iVV„^/^(r - N) + r^(3r^ + 24N' + f) 

+5A^r2(3r2 + f) + rN^{-27N'^ + 7/^) + ?>N'^{-hN^ + f)]. (63) 

It is worth to note that the Taub-NUT solution has not curvature singularity at r = A^ for 
B = CP , while for B = S^ x S^ the metric has curvature singularity. Since Fn{r) depends 
on Vn, the fourth condition for NUT solution restricts the allowed values oi Vn- As we 
discussed in Sec. 11111 one may obtained the critical value of Krit by solving the system of 
two equations ^I7\. To be more clear, we obtain Vcrit for asymptotically flat (/ -^ oo) NUT 
solutions. The system of two equations ^T7\ becomes 

(r + 3N){r + N)^ - 12N\r - N)V^^^, = 0, 
(3r + 7N){r + N) - 12N^V^^-^, = 

with the following solution for V^rit 

^„,..^«^^.l,, (64) 

D 

For arbitrary values of /, one may find the critical value of V numerically. For / = 1 and 
A^ = 1 the critical value of potential which is obtained by the above method is Vcrit = 4.67. 
This can be seen in Fig. |3] which shows the function Fn{r) as a function of r for various 
values of V including V = V^rit- The electric charge and potential are 

Qn = I27iq, (65) 

^" - ^- (66) 

The Euclidean actions can be calculated by use of Eqs. (0), dHl) and Q as 

if-AN^-m^V^), (67) 



3/2 
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FIG. 4: Fnir) versus r for F = 3 < Fcrit (point), V = Fcrit = 4.67 (line) and F = 6 > Fcrit 
(bold-line). 



It is notable to mention that the electromagnetic field will cause a power law divergence 
IGiiV'^N'^r in action which is removed by the counterterm (j^I). The stress-energy tensor 
corresponding to the counterterm (jH)) in six dimension is 

rp{em) _ ^ { rpc rp ^ rp r?cd^ \ 

J-ah - ^K.i' a^ cb - ^r cd^ gab) 

The mass due to this counterterm is 



(68) 



M^r^ = levriVW, 



(69) 



which removes the power law divergence of the mass due to electromagnetic field. The total 
mass is 

Mn = ^-nmn (70) 

The Smarr-type formula for mass is 



Mr 



967rZ3/2/2 



(71) 



where Z = N"^ is a function of Q and S which may be obtained by the following equation 



16 ^ 



The entropy and specific heat are 



S^ = 2^lJ:L^20N^-3l''-5AfV^), 
a^^^(/^-10iV^-9/V^). 



(72) 

(73) 
(74) 
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As one can see from Eq. UM for \V\ < (6/)-^/(4073)iV^^2/2 and \N\ > VoJEl, the entropy 
is positive, but C„ is negative in this range. Thus, S'„ and C„ are not positive simultaneously, 
and the NUT black hole is completely unstable in this ensemble as in the case of uncharged 

n 

solution )30l ]. 

B. Bolt Solutions 

The parameters nib, Qb and Pb are 

mb = -^[Sr^b + (^' - l^Ny^ - 3N\2f - 15Ny? - 3N\P - 5N^) (75) 

3rrb 

-^(^t + 3^M + 15iVV,2 - SN% (76) 

Vjrt - QN'rl - 3N') 
g. = ' (77) 

^'^ " 2(5r2 - 5m + P)rl + 9(r2 - m)PV^- ^^^^ 

where Vb is real root of the following equation 

30N4 - 2frl + 3(2/2 _ ^q^2 ^ gpy2^j^^2 _ 27Ar3/2y2 _ Q_ ^7g^ 

Equation ()54|) has two real solutions greater than N provided A^ < A^max and one real 
solution Tb > N for A^ = A'max where A^tnax is the smallest root of the following equation 

9 X (10)^A^^° + (10)^ X (324^2 - 72)/2A^^ + (10)^ x (43741/^ - 648^^ + 217)/^ A^^ 
+ (10)=^ X (26244^6 + 5832^^ - 999^^ - 294)/6A^^ + (10) x (59049V* + 52488^^ 
+14823 V^ + 2052^2 + 156)/* A^^ - (729V^ + 486V^ + 811/^ + 8)/^° = 0. (80) 

As in the case of 4-dimensional solutions, A'max decreases as V increases. 
The electric potential is 

^<> = - rl - 6N^rl - 3^' ''^^^ 

The finite action, entropy and heat capacity can be obtained as 

_ 87r=^ rg[-3rg + {P + hN^)4 + 3N''{hN^ - 2f)rl + 3N\hN^ - f)] 
^ ~ 3 2(5r2 - 5Ar2 + /2)^2 ^ 9(^2 _ n2^12v2 

I 1 o,2,2..2 {rl - 15ArV,^ - 21N^tI + 3N^) 

2(5r2 - 5Ar2 + /2)^2 + 9(^2 _ f^2^^l2y2 ' ^^^) 
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FIG. 5: Entropy (bold line) and heat capacity (solid line) of 6-diniensional bolt solution for / = 1 
and y = in the allowed range of N. 



s, 






2;2t^2 



+367r'rF 



2r2(5r2 - 5N^ + P) + 9(r2 - N^)PV^ 



(83) 



2„4 



+220N^4 - lOl^lN^ - bONrf + 3rlf]}-^[-3600N^H^rt - AAON^rfl^ + gSOOiV^/V^ 
+90QN^^rlf + 20Nrl'^l^ - IbQN^ll'^ - USSON'^frf - 288iVVf /^ + 18000A^Vf /^ 
+1152A^V?/^ - 1680A^^/VP - AGOONH'^rl^ + 160A^VH/^ - 8280N\fP - breN^rll"^ 



+324N^l^rt + 840N^r^l^ + m20N^rlf + 200m^r}'^ + ISON^IH^ - GOrlH^ - UrlH* 



+23100N^rl + 500Nrl^ - 204800A^\^° + 62500N^rl^ + iOON^rl^ + 9900A^^V^ 
+6000A^^V^ + (31536A^V^/^ - 20A12NH^r^ + 7938N^rll'^ + dOONl^l'^ - QlbGNVj"^ 
-12312N^r^J* - 626AN^rl^l^ + 378N^\ll^ + 53QAN^rll^ + 180A^rf /^ - 19AAN^^rll^ 
+113A0N^\ll^ - 386100iV^r^/2 + 17A6N\IH^ + 159840iV^/V^ + 116640A^^^/V^ 
-52380N^rl°l'^ - 3A5Q0NH^rl'^ - 5ArlH^)V^ + (19440A^V^/^ - 10Q92N^rl°l'^ 
-14823ArV^/^ - 117369Ar^/V + 37908N^^rll'^ + AObNrlH^ + 2187N^H'^)V% (84) 

Figures El and m show the entropy and specific heat as a function of A^ for / = 1, V^ = 0, 
V = 1 and l^ = 3 up to A^max for outer horizons of the boh solutions respectively. As one 
can see from these figures, the bolt solution is stable in the range < A^ < A^st < ^max, 
while A'st is less than A'max for V ^ 0, and equal to A^max for V = 0. That is, the uncharged 



19 



1i!+06- 
80(1000 z 
60(1000 
40(1000 
20(1000 z 





-20(f000 
-40(i000 




0.01 0.015 0.02 ,/fl.025 0.03 



FIG. 6: Entropy (bold line) and heat capacity (solid line) of 6-dimensional bolt solution for / = 1 
and y = 1 in the allowed range of N. 




FIG. 7: Entropy (bold line) and heat capacity (solid line) of 4-dimensional bolt solution for I = 1 
and y = 3 in the allowed range of N. 

solution is stable in the whole allowed range < A^ < A^max? while the stable phase becomes 
smaller as V increases. 



VI. THERMODYNAMICS OF {2k + 2)-DIMENSIONAL CHARGED TAUB-NUT 
SOLUTIONS 



Using Eqs. Q, (jHI) and ©, the Euclidean actions can be calculated as 
/^^ = -J(fc)[/2 _ 2kN^ + 2{-l)''+\2k - 1)H^V% 
where J{k) is a function of k as 

J{k) 



iA7r)'+\k + im-k)Tik + ll^,, 



167r3/2/2 



Also the total mass may be calculated as 



\k-l 



Mn = t: mn. 



(85) 



(86) 



(87) 
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As in the case of uncharged solutions [30| , the total angular momentum is zero. The entropy 
can also be obtained through the use of Gibbs-Duhem relation (Q), where Ci = Q and Fj = $ 

as 

S^ = -J{k) [{2k - 1)P - 2k{2k + 1)N^ + 2{2k - ifl^V^] (88) 

The Smarr-type formula for mass versus extensive quantities S and Q may be written as 

M^ = M4>r)'-'B«:)Z<-'^'P - 2(k + 1)Z| + (^'^I'U, l r(t + 1)^! - l) '''?'' '»"' 
where Z is the solution of the following equation: 

UTxf+\k+l)Ti}:-k)Tik + l) ^, 

Sn + 16^3/,^, -m - 1)1' - 2k{2k + l)Z 

,(2fc-i)/^ r(fc-|) ^.^^ , 



TmYQ'}z' = 0. (90) 



27r '(47r)'=-ir(A; + l)Z(fc-i/2)J 

One may then regard the parameters S and Q as a set of extensive parameters for the 
mass M.{S,Q) and define the intensive parameters conjugate to S and Q, which are the 
temperature and the electric potential respectively. One obtains 

'dMn\ 1 

Q 



^" ' dS Jr> 4:{k + l)7rN' 



Equations ()9H1 show that the the quantities T„ and $„ coincide with those which was 
calculated in Sec. pil|) . Thus the thermodynamic quantities calculated in Sec. pil|) for 
Nut solution, satisfy the first law of thermodynamics. The specific heat capacity may be 
obtained as 

Cn = 2J{k)[k{2k - l)f - 2k{k + l){2k + 1)N^ - 2{k - l)(2k - lffV\ (92) 

Since J{k) will produce a minus sign for odd fc, Eqs. ()88|) and (|U^ show that the entropy 
and heat capacity are positive for odd k provided 

2fc(2fc + l) ,,2 ..., ..2,2..2 .,2 ^ 2fc(fc + l)(2A: + l)iV^ + 2(A:-l)(2fc-l)3/^V-^ 
{2k -I) ^^ ~ ^^^^ ~^^^^ < ^ < A;(2A;-1) 

Hence, for odd fc > 1, the NUT solutions will be thermally stable in the above range, which 
becomes increasingly narrow with increasing dimensionality and wide with increasing V . 
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For k = 1, as we see in Sec. IIV[ there exists another stable phase. For even k, however, no 
minus sign is produced by J{k), meaning that the entropy will be positive for 

[jdrh J- j 

and the specific heat will be positive for 

2 2k{k + l){2k + 1)N^ + 2{k - l){2k - ifl^V^ 
^ k{2k-l) 

Since the second value will always be larger than the first, this means there is no range in 
which both will be positive, and thus for all even k, the NUT solutions will be thermally 
unstable. 

VII. CONCLUSION 

In this paper, we first, constructed asymptotically AdS Taub-NUT/bolt solutions in 
{2k + 2)-dimensional Einstein-Maxwell gravity with curved base spaces. We found that 
the function F{r) of the metric does not depend on the specific form of the base factors on 
which one constructs the circle fibration. In the presence of electromagnetic field, there exist 
two extra parameters, in addition to the mass and the NUT charge, namely; the electric 
charge q and the potential at infinity V. We found that in order to have NUT charged 
black holes in Einstein-Maxwell gravity, in addition to the two conditions of uncharged 
NUT solutions, there exists two other conditions. The first extra condition comes from the 
regularity of vector potential at r = A^ which gives a relation between q and V. Indeed, 
the existence of the parameter V enables us to get a regularity condition on the one-form 
potential which is identical to that required to obtain a NUT solution. If one of these 
parameters vanishes then the other one should be equal to zero and the solution reduces to 
the uncharged solution. The second extra condition comes from the fact that the horizon 
at r = A^ should be the outer horizon of the black hole. Indeed, Einstein-Maxwell gravity 
admits NUT black holes provided the potential parameter is less than a critical value V^rit? 
which may be obtained by solving the system of two equations ()27|). In any dimension larger 
than four, the mass parameter m which is fixed by these four NUT conditions depends on 
the fundamental constant / and parameters A^ and V (or g), while in 4-dimensions the mass 
parameter does not depend on V (or q). We also found that the NUT solutions have no 
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curvature singularity at r = A^ when the metric of the base space is chosen to be CP , 
and have curvature singularity for other curved base spaces. Then, we obtained the bolt 
solutions of Einstein-Maxwell gravity in various dimensions, and gave the equations which 
can be solved for the horizon radius of the bolt solution. We found that in order to have 
bolt solution, the NUT parameter should be smaller than a maximum value A^max which 
decreases as V increases. 

Second, we investigated the thermodynamics of NUT solutions in any dimensions. We 
calculated the conserved quantities and the Euclidean actions of the NUT/bolt solutions 
through the use of counterterms renormalization procedure. Also we obtained the charge 
and electric potentials of the solutions in an arbitrary dimension, and calculated the 
entropy through the use of Gibbs-Duhem relation. We obtained a Smarr-type formula for 
the mass as a function of the extensive parameters S and Q, calculated the temperature 
and electric potential, and showed that these quantities satisfy the first law of thermody- 
namics. Then, we studied the phase behavior of the charged NUT solutions in canonical 
ensemble by calculating the heat capacity, and found that the NUT solutions are not 
thermally stable for even k, while there exists a stable phase for odd k, which becomes 
increasingly narrow with increasing dimensionality and wide with increasing V. We also 
studied the phase behavior of bolt solutions in 4 and 6 dimensions in canonical ensem- 
ble and found that these solutions have a stable phase, which becomes smaller as V increases. 
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